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Abstract 

A detailed numerical stability analysis of the static, spherically symmetric globally reg- 
ular solutions of the Einstein- Yang-Mills equations with a positive cosmological constant, 
A, is carried out. It is found that the number of unstable modes in the even parity sector 
is n for solutions with n = 1,2 nodes as A varies. The solution with n = 3 nodes exhibits 
a rather surprising behaviour in that the number of its unstable modes jumps from 3 to 
1 as A crosses (from below) a critical value. In particular the topologically 3-sphere type 
solution with n = 3 nodes has only a single unstable mode. 
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As we have learned during the last few years, the coupling of nonabelian gauge fields 
to Einstein's gravity leads to a rather large class of solutions of potential physical interest, 
such as "gravitating glueballs" of Bartnik-McKinnon 0, and "hairy black holes" |21 El H|- 
These globally regular solutions are static, spherically symmetric, and they form a discrete 
family labelled by the number of nodes, n, of the (single and purely "magnetic") Yang-Mills 
amplitude. Furthermore the above solutions have been generalized to a continuous family by 
allowing for a positive cosmological constant, A (Sj. These "cosmological" solutions exist only 
for a limited range of the cosmological constant, < A < A reg (n), n being the node number of 
the Yang-Mills amplitude. 

The important problem of stability has been studied mostly in linear perturbation theory, 
with the conclusion that all of these solutions are unstable El El 03 El- The spectrum 
of small fluctuations can be divided with respect to parity into even ("gravitational") and 
odd ( "sphaleronic" ) sectors. Within the most general spherically symmetric ansatz, mostly 
numerical results give convincing evidence that a gravitating glueball or black hole with n 
nodes possesses exactly n unstable modes in the even parity sector 00, and in Refs. 015] 
analytic arguments have been presented that there are also precisely n unstable modes in the 
odd parity sector. 

The aim of this paper is to point out that in the even parity sector the number of unstable 
modes of the cosmological solutions shows an unexpected dependence on the cosmological 
constant, A, contrary to a previous claim P|. Since the geometry of the cosmological solutions 
changes substantially for a certain value of A = A cr i t (n), as then an equator develops, i.e. the 
surface area, 4irr 2 , of the invariant 2-spheres becomes stationary for a certain value of r and 
then it starts to decrease. A careful numerical analysis reveals that the number of unstable 
modes in the even parity sector is 1 resp. 2 for solutions with n = 1 resp. n = 2 nodes 
independently of A as it varies from to A reg (n). The number of unstable modes of the n = 3 
nodes solution exhibits, however, a rather surprising behaviour. As A tends to a certain value, 
Amerg ~ 0.2658, two eigenvalues of the unstable modes approach each other, actually merging 
at A merg and then they simply disappear from the spectrum. In particular the globally regular, 
spatially compact and topologically 3-sphere-type solution (where A = A reg (n)) with n = 3 
zeros has only a single unstable mode. This is a somewhat unexpected phenomenon, since the 
linearized perturbation equations can be transformed to a standard 1 dimensional Schrddinger 
eigenvalue problem for the perturbation, SW, of the Yang-Mills background for a particular 
choice of the coordinates (Schwarzschild gauge). The unstable modes are then nothing but 
bound states of the corresponding Schrddinger equation. Now for such a Schrddinger problem 
one does not expect that the values of bound-state energy levels could even coincide as a 
single parameter (A) is varied, let alone disappear from the spectrum without crossing the zero 
energy threshold. The Schwarzschild gauge is, however, singular when an equator is present, 
which is always the case for A £ [A cr it (n) > A reg (n)] . In that case the potential in the associated 
Schrddinger equation becomes so singular that the standard theory of self adjoint operators is 
no longer applicable. It is quite remarkable that the number of unstable modes can actually 
change as the cosmological constant varies and it would be clearly desirable to achieve a better 
understanding of this phenomenon. 

The most general spherically symmetric line element can be written as 

ds 2 = e 2 ^dt 2 - e 2X ^dR 2 - r 2 (R, t)dQ 2 , (1) 
and the minimal spherically symmetric Ansatz for the YM field is given by 

A = (T 1 W(R, t)d6 + T 2 W{R, t) sin Qdsj> + T 3 cos Od<fy , (2) 

where the T l {i = 1,2,3) denote an orthonormal set of Lie algebra generators of SU(2). For 
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the above Ansatze the (suitably rescaled) Einstein equations take the form: 

1-U-2V- e" 2 V(r 2 )' + e- v ({r 2 ) it e- v ) tt = , (3a) 
1 + U - 2V - e- x (e- x (r 2 )'y + e- 2u \ t {r 2 ), t = , (3b) 

e- u {r, t e-»), - e- A (r'e" A )' - e^- x [(e^ x ru'Y - (e^rX,),] - ^- = , (3c) 

r ' +r 'A ( - r / + S! = , (3d) 
r 

together with the Yang-Mills field equation for the amplitude W: 

{e u - x W'y - (e x - u W :t ), t - e u+x W ^ 2 ~^ = , (4) 

where 2 

V = {l ~ W 2 ] + , U = e- 2 »{r 2 - 2W 2 ) + e~ 2X (r' 2 - 2W> 2 ) , 
2r z 2 

and r jt := jfe, while prime denotes 

Following Ref. the static field eqs. are conveniently written as 

r = N, (5a) 
N W 2 

N=( K -N)—-2—, (5b) 

k = \l - k 2 + 2P^ 2 - 2Ar 2 j fr , (5c) 

W 2 - I W 
W = W-— (k-N)— , (5d) 

where from now on / := 4^ := e~ x f, and JV, k are defined as N := f, k := N + ri>. The (i?i?) 
component of the Einstein eqs. (jHK) yields the constraint : 

1 + iV 2 + 2W 2 - 2V = 2kN . (6) 

Solutions of Eqs. © with a regular origin have the following power series expansion 

r = <7-i(2& 2 + fy 3 + 0(a 5 ), k = 1 + (26 2 - A ) ff 2 + 0(a 4 ) , i/ = (26 2 - £ )a + 0(<r 3 ) , 

iV = 1 - (26 2 + |)a 2 + 0(<r 4 ) , = 1 - bo 2 + 0(a 4 ) , 

(7) 

where 6 is a free parameter. We consider solutions with a cosmological horizon, characterized in 
our notations by the vanishing of N(<Jh) and a simple pole of k(<t) at a = a^. The corresponding 
series expansions of the functions near the cosmological horizon are given as: 



r = r h + r[ 2) x 2 + 0(x 4 ) , k = — + k^x + 0(x 3 ) , u = v h + ln(x) + 0(x 2 ) 

N = N^x + C»(x 3 ) , W = W h + W ( h 2) x 2 + 0(x 4 ) , 
with x = ah — a, and 



(8) 



3 V r h 2 



h 2r h \ r 2 h I h 4 r 2 



(9) 



where Wh, Vh an d r~h are free parameters. 

As found in Ref. [Jj the cosmological solutions naturally fall into three qualitatively 
different classes as the cosmological constant varies from A = to A = A reg (re). All of these 
solutions may be indexed by the number of zeros, n, of the Yang-Mills amplitude, W. For 
A = the space-time corresponding to the globally regular Bartnik-McKinnon solutions is 
asymptotically flat. As soon as A > the globally regular cosmological solutions correspond to 
asymptotically de Sitter space-times (with a cosmological horizon), as long as the cosmological 
constant does not exceed a critical value, i.e. < A < A cr it(n), constituting the first class. 
The geometry of the solutions changes when A exceeds A CT [t(n), since then an equator develops 
outside of the cosmological horizon, i.e. the radius of the invariant 2-spheres, r(cr), reaches a 
maximum and then it decreases down to zero where a geometrical singularity develops. The 
equator is located outside of the cosmological horizon for A < A*(n). As A further increases, 
A*(n) < A < A rcg (n), the position of the equator moves inside the horizon. Such solutions 
with an equator correspond to a singular geometry of the "bag of gold" type. Finally for a 
special value A = A reg (n) the position of the horizon coincides with the singularity and then 
remarkably, there exists a globally regular solution corresponding to a space-time whose spatial 
sections are compact. No regular solution seems to exist outside this range of the cosmological 
constant. For n = 1, this solution is analytically known 



K = N = W = cos(a/V2) , r = V2sm(a/V2) , A reg (l) = 3/4, (10) 

corresponding to a static Einstein universe whose spatial section are 3-spheres, with a constant 
energy density due to the Yang-Mills fields. 

To study the linear stability problem of the above solutions, one introduces time depen- 
dent perturbations of the form 

f(a,t) = f(a) + 6f(a)e i » t , (11) 

where now / stands for any of the functions {r, N, k, W}, and one linearizes the field eqs. (j2K-d) 
around a static background, f(cr). In this way one obtains the following set of linear differential 
eqs. for the perturbations: 

$ r = 1^—5r- ^Ksw + NSX, (12a) 

K — IV W f 3k \ 

5k=(k- e~ 2v u?r)5X - 2- 5k + A—5W - 4A + — + e~ 2i V 5r , (12b) 

r r \ r / 

k — N ■ ( 3W 2 - 1 2W 2 9 „ 9 \ WW 
5W = 5W + - ^ - e~ 2u uj 2 5W - 2—5r 5k 

/ K — N ■ \ 

+ 12W + W 5\ + W5X , (12c) 



r 

with the linearized constraint (pji): 

N5k -(k- e - 2u to 2 r)5r + 2W5W - 2W5W + 2W 2 5\ = . (13) 

We remark that Eq. (|12b ) is nothing but the (t, R) component of the linearized Einstein equa- 
tions. The system ()12|) is invariant with respect to the following gauge (diffeomorphism) trans- 
formations: 

5X — ► 5X + F , 5k — ► 5k + (k - e~ 2v uj 2 r)F , 

(14) 

5r — > 5r + NF , 5W — > 5W + WF , 

where F is a function of the variable a. As already mentioned, due to the presence of an 
equator the stability analysis of the cosmological solutions is more complicated than that of 
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asymptotically flat ones. When there is no equator there exists a particular regular gauge - the 
Schwarzschild gauge - defined by 5r = 0, in which the gravitational degrees of freedom can be 
eliminated from the equation of the Yang-Mills perturbation amplitude in Eq. ()12b ) . Indeed, 
in this gauge one obtains just a standard Schrodinger equation: 



? (fbW 



+ V S 5W = uj 2 5W , (15) 



where the potential, given by 

V s = <?" "" = - +Ae»\e v ^-\ , (1.0) 




becomes singular when the solution has an equator, V s oc 1/(<t — o" eq ) , due to the vanishing of 



N. The unstable mode of the analytically known 3-sphere background (|10|) corresponds to the 
following solution of the Schrodinger equation (|15j): 

5W = sn(a/V2) tan(<x/\/2), u? = -1. (17) 

Note that the solution l|17jl is not even square integrable because of its divergence at the 
equator a eq = ir/y/2. According to the classical theory of self-adjoint operators one should 
impose boundary conditions at the equator which would, however, just exclude the solution 
Q17|) (see Ref. 9 for a discussion on this point), so it seems that the correct interpretation of 
the singular Schrodinger equation ()15|) in the present context is beyond the standard theory. It 
is clear for physical reasons that one cannot impose boundary conditions for the perturbations 
at the equator, since the singularity there is an artifact due to the gauge choice. 

On the other hand one can simply choose a globally regular gauge where the equator is 
a regular point of Eqs. (|12|1. The only disadvantage of such regular gauges is that one does not 
obtain a Schrodinger- type perturbation equation, but a somewhat unusual eigenvalue problem. 
It is then not easy to interpret the eigenfunctions and give a meaning to the number of nodes 
as in the case of a Schrodinger equation. Nevertheless we find it both conceptually clearer and 
also better suited for our numerical procedure to work in a globally regular gauge, which we 
have chosen to be 5X = 0. 

The solutions of Eqs. (|12[) corresponding to locally regular perturbations admit the fol- 
lowing power series expansions at the origin: 

Sr(a) = ^-5W 2 a 3 + 0(cr 5 ) , 5k{cj) = -4b5W 2 a 2 + 0(o A ) , 5W(a) = 5W 2 a 2 + 0(a 4 ) , (18) 
and near the horizon 

5r(x) = 5r h x a (x 2 +0(x 4 )) , Sk(x) = 5K h x a {x+0{x 3 )) , 8W{x) = 5W h x a {l+0{x 2 )) , (19) 
with the following definitions 

a 2 = - e - 2 ^ 2 , 5r h = -^^-5W h , 6 Kh = -^~-^5r h , (20) 

r h (a + l) (a + 3) 

where 5W 2 , u 2 and SWh are free parameters. For the background (|TU|) the unstable mode can 
be written in the 5X = gauge as: 

8W = -5k = (cr/\/2) sm(a/V2) , 5r = ^/2 sin(CT/V2) - a cos(ct/\/2) . (21) 

We have numerically integrated Eqs. (|12[) by a fifth order adaptive step-size Runge-Kutta 
method from both the origin and the horizon subject to the boundary conditions (|18|) and (|19|) 
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in a given static background solution, and by varying the free parameters we have matched the 
solution at an intermediate point (shooting to a fitting point). 

Our main numerical results are presented on Figures 1 (n = 1, 2, 3) and in Tables 1-4. The 
A-dependence of the unstable mode of the n = 1 background solution displayed on Fig. l(n = 1) 
agrees with that of Ref. The main new behaviour of the unstable modes is exhibited on Fig. 
1 (n = 3): the largest eigenvalue unstable mode with two nodes of the n = 3 node background 
solution - corresponding to the highest energy bound state wave function in the Schrodinger 
description - looses one of its nodes at A merg ~ 0.2657847 where its eigenvalue coincides with 
that of the single node bound state. At A merg the coinciding eigenvalues are w| = w\ « —69. 
For A > A merg these modes with nodes completely disappear from the spectrum, so only a 
single unstable mode without nodes (the lowest lying bound state) survives. It is worth noting 
that the lowest lying eigenvalues of the bound states for the n = 2 and in particular for the 
n = 3 node background solution exhibit quite large variations in a rather small A interval close 
to A = 0. For example the energy of the lowest lying bound state in the n = 3 background 
drops from ~ —3000 to ~ —190000 as A varies from to ~ 0.02. For backgrounds with n > 3 
this phenomenon is getting more and more pronounced, making numerical integration quite 
difficult. 

Preliminary investigations |13j suggest that in fact for all n > 3 node solutions there is 
such a A merg (n) where precisely two negative eigenvalues collide, leading to a diminution of the 
unstable modes by two. Therefore it is natural to conjecture that the number of instabilities 
of the n > 3 node topologically 3-sphere-type solution in the even parity sector is in fact n — 2. 

Finally for completeness we have also computed the eigenvalues of the unstable modes 
in the odd parity sector, in a regular gauge for solutions with A = and for the topologically 
3-sphere solutions. Omitting all computational details here, we just state that our results show 
that the number of unstable modes for both A = and A = A rcg (n) is n (up to n = 3), in 
agreement with the analytical arguments of Ref. [Hj. The numerical eigenvalues are provided 
in Table 4. 
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Table 1: Numerical values of the parameters for some n = 1 solutions. 



A 


b 




w h 


^0 





0.45371627 


oo 


-1 


-3.284026 


0.001 


0.45358471 


53.924708 


-0.983515 


-3.283974 


0.01 


0.45234499 


16.431260 


-0.947872 


-3.287678 


0.1 


0.43582245 


4.441700 


-0.859362 


-3.162423 


0.5 


0.32415673 


0.956828 


-0.904637 


-1.649222 


0.7 


0.26406135 


0.208989 


-0.990741 


-1.097542 


0.75 


0.25 





-1 


-1 



Table 2: Numerical values of the parameters for some n = 2 solutions. 



A 


b 




w h 


-i 2 


^0 









0.65172553 


oo 


1 




-17.7898 




-93.9115 






0.001 


0.65357281 


53.775678 


0.837449 




-20.0993 




100.6350 






0.01 


0.66188220 


16.271892 


0.560250 




-39.3607 




153.0908 






0.1 


0.63599416 


4.224264 


0.332182 




-76.9804 




222.8341 






0.3 


0.49283156 


1.244954 


0.647999 




-13.2317 




-17.9967 






0.35 


0.44326287 


0.419923 


0.938400 




-5.7734 




-9.5080 






0.36424423 


0.42959769 





1 




-4.7046 




-8.0863 






Table 3: Numerical values of the parameters for some n = 3 solutions. 






A 


b 


rh 


w h 




W jf 


^0 





0.69704005 


oo 


-1 




-53.3132 




-390.5425 


-2946.8374 


0.001 


0.70174503 


53.754661 


-0.321763 




-370.5305 


-1235.9470 


-14402.84 


0.01 


0.70073671 


16.257493 


-0.119067 




-764.4160 


-9097.367 


-94777.96 


0.1 


0.66137928 


4.206730 


-0.075888 




-493.7848 


-9938.45 


-89162.65 


0.26578472 


0.54488550 


1.351120 


-0.444434 




-69.0027 




-69.0678 


-153.7337 


0.29 


0.51312783 


0.347933 


-0.947303 












-35.7286 


0.29321764 


0.50882906 





-1 












-30.7787 
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Table 4: Eigenvalues of the unstable modes in the odd parity sector. 



n 


A = 


A = Areg(n) 


1 


-3.872828 


-1 


2 


-24.083467 


-3.665216 




-82.480704 


-4.619788 


3 


-73.653692 


-7.938386 




-325.34471 


-9.198014 




-3005.4050 


-33.599200 
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